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ABSTRACT 

In this paper the synthetic period- luminosity (P-L) relations in Spitzer's IRAC bands, based on a 
series of theoretical pulsation models with varying metal and helium abundance, were investigated. 
Selected sets of these synthetic P-L relations were compared to the empirical IRAC band P-L relations 
recently determined from Galactic and Magellanic Clouds Cepheids. For the Galactic case, synthetic 
P-L relations from model sets with (Y = 0.26, Z = 0.01), (Y = 0.26, Z = 0.02) and (Y = 0.28, 
Z = 0.02) agree with the empirical Galactic P-L relations derived from the Hubble Space Telescope 
parallaxes. For Magellanic Cloud Cepheids, the synthetic P-L relations from model sets with {Y = 
0.25, Z = 0.008) agree with both of the empirical Large Magellanic Cloud (LMC) and Small Magellanic 
Cloud (SMC) P-L relations. Analysis of the synthetic P-L relations from all model sets suggested that 
the IRAC band P-L relations may not be independent of metallicity, as the P-L slopes and intercepts 
could be affected by the metallicity and/or helium abundance. We also derive the synthetic period- 
color (P-C) relations in the IRAC bands. Non- vanishing synthetic P-C relations were found for certain 
combinations of IRAC band filters and metallicity. However, the synthetic P-C relations disagreed 
with the [3.6] — [8.0] P-C relation recently found for the Galactic Cepheids. The synthetic [3.6] — [4.5] 
P-C slope from (Y = 0.25, Z = 0.008) model set, on the other hand, is in excellent agreement to the 
empirical LMC P-C counterpart, if a period range of 1.0 < log(P) < 1.8 is adopted. 
Subject headings: distance scale — stars: variables: Cepheids 



1. INTRODUCTION 

The mid- infrared Cepheid period-luminosity (P-L, also 
known as the Leavitt Law) relation will be important in 
the James Webb Space Telescope (JWST) era, as it holds 
the promi se of deriv ing the Hubble const ant at the ^ 2% 
level (Frce dman fc Madorc 2010; Frcedm aiTet al.ll2011[ ). 
Motivated by this, the Spitzer's IRAC band (3.6/im, 
4.5/im, 5.8yLtm and 8.0/^m) P-L relations were de rived for 
Cepheids in our Galaxy (iMarengo et al.ll2010D. i n the 
Large Magellanic Cloud f LMC. iFre edman et all 120081: 
Ngcow & Kanbur 2008; Mad ore et al . 2009; Ngcow et ail 
20091: IScowcroft et al. 2011), and in the Small Magellanic 
Cloud (SMC; |Ngeow & Kanbur 2010). 

The slopes of the IRAC band P-L r elations are ex- 
pected to be insensitive to metallicitv (IFreedman et al.l 
l2008l:lRiedman fc MadordlMoHlYeedman et al.ll201lD . 
The empirical slopes derived from a small num- 
ber of Gala ctic Cepheids that possess parallax dis- 
tances (|Marengo et al.ll2010f) . and the Magellanic Cloud 



Cepheids based on the OGLE -HI data (see lNgeow et al.l 
I2OOI INgeow fc KanbiullMo L for more details) are all 
consistent with each other. However, these slopes do 
not agree with the slopes derived from Galactic Cepheids 
that are based on the infrared surface brightness (IRSB) 
method ([Marengo et al.J | 2010f) and fro m a smaller num- 
ber of LMC'Cepheids'^^adore^^ Therefore, 
the aim of this paper is to compare these empirical P- 
L relations to the synthetic IRAC band P-L relations 
based on a series of theoretical pulsating models at vari- 
ous metallicities, and investigate the sensitivity of IRAC 
band PL relations to metallicity. 

Brief description of the pulsation models is given in 
the next section, and the synthetic P-L relations based 
on these models are presented in Section 3. We com- 
pared these synthetic P-L relations to their empirical 
counterparts and investigated the possible metallicity de- 
pendence of the synthetic IRAC band P-L relations in 
Sections 4 and 5, respectively. In Section 6, synthetic 
period-color (P-C) relations in the IRAC bands were also 
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derived. Discussion and conclusion are presented in Sec- 
tion 7. 

2. THE PULSATIONAL MODELS 

The synthetic P-L relations adopted in this paper are 
based on extensive and detailed sets of nonlinear, pulsa- 
tion models including a non-local, time-dependent treat- 
ment of the coupling between pulsation and convection. 
These models allow us to predict not only the peri- 
ods and the blue boundary of the instability strip, but 
also the pulsation amplitudes, the detailed light and ra- 
dial velocity curve morphology, and the com plete topol 
ogy of the strip, including the red edge ([Bono et aU 



Y=0.25 2=0.008 



[Tool iFiorentino et al.|[200a iMarconi et al.| [2005. 201 

For each chemical composition and mass, an evolution- 
ary mass-luminosity (M-L) relation was adopted (see 
IMarconi et al.l [20051 for details) and a wide range of ef- 
fective temperature was explored. We note that even if 
the effect of varying the M-L relations has been investi- 
gated for specific chemical compositions (see, for exam- 
ples, [Bono ct al. 1999, 2000; Caputo et al. 2005) in this 
analysis we assume for all the chemical composition a 
canonical M-L relation, neglecting both mass loss and 
overshooting during the previous H burning phase. This 
is a limitation of the adopted model sets as the above 
mentioned phenomena affect the P-L relations and the 
corresponding distance determinations (see Section 7). 
In total, 17 sets of model with varying helium {Y) and 
metal {Z) abundance are considered in this 



of them are the same model sets presented in lBono et al.l 
(HoiO). 

From the resulting theoretical instability strips and 
the relations connecting the periods to the intrinsic stel- 
lar parameters, synthetic P-L relations have been con- 
structed. To this purpose, we populated the predicted 
ins tability strip by adopt ing the procedure suggested 
bv iKennicutt et all (|1998[ ). In particular, 1000 pul- 
sators were uniformly distributed from the blue to the 
red boundary of the instability strip, with a mass law as 
give n by dn/dm = over the mass range 5 - llAf© 
(see lCaputo et al.l l2000L for further details). In order to 
translate the pulsational properties of the investigated 
Cepheid models in the Spitzer IRAC bands, we have 
directly convolved the predicted bolometric light curves 
with the Spitzer filter profi les using the genera l integral 
equation (see, for example. iGirardi et al.il200"2[ ): 



-2.5 log 



(1) 



where Sx is the IRAC Spectral Response Curv£[3, fx is 
the stellar fiux (that corresponds to model atmospheres 
of known {Tcs, [AI/H],\ogg)), /° is the model spectrum 
of Vega. Concerning the model atmospheres, we have 
adopted the homogeneous set of updated ATLAS9 Ku- 
rucz model atmospheres and synthetic fluxes (new-ODF 
models j3- 

3. THE SYNTHETIC IRAC BAND P-L RELATIONS 



^ Available at |htt p: //ssc ■ spitzer ■ caltech. edu/irac/| 
^ Available at http://kurucz.harvard.edu/grids.htnil 
http : //wwwuser . oat . ts . astro . it/castelli/grids .html 




Fig. 1. — Example of the synthetic IRAC band P-L relations, 
constructed from models with Y = 0.25 and Z = 0.008. 



Figure [T] shows an example of the synthetic IRAC band 
P-L relations from one of the model sets described in the 
previous section. All the synthetic P-L relations were fit- 
ted with the form of Mirac = a+bxlog{P), restricted for 
pulsat ors within the pe riod range of 0.4 < log(P) < 2.0 
(as in lBono et aIll2010D . The fitted P-L slopes (b) and 
intercepts (a) for each of the models were summarized 
in Table [T] and [U respectively. A majority of the model 
sets do not have pulsators with log(P) < 0.4. For a few 
of the model sets which have a small number of short- 
period pulsators, the differences of the P-L slopes and 
intercepts between the P-L relations derived from full 
period range and those given in Tables [1] and [2| do not 
exceed 0.01. In Figure[2l the slopes of the synthetic P-L 
relations in various bands were compared for six of the 
pulsating model sets (or chemical compositions), where 
t he linear BVI JK P -L slopes were adopted from Table 2 
of lBono et al] (|2010lFI . As expected, the slopes of the P-L 
relation mon otonically d ecrease fr om B to K ban d (see, 
for ex ample. iMadore &: Frccdma nlll991l: [B crdniko v et"ail 
11996"; 'Caputo et al.' '2000; Fiorentino ct al. 2002|, 120071 : 
Freedman ct al. 2008; Ngeow & Kanbur 2008), and "flat- 
ten out" in the mid-infrared. However, the 4.5^m P- 
L relations show a slight increase in their slopes when 
compared to the "flatter" slopes defined from 3.6/xm and 
8.0/im band P-L slopes. This slight increase of the 4.5/im 
P-L slopes, and in some extent to the 5.8/im band P-L 
slopes, may be explained due to the presence of CO ab- 
sorption features shown in the ~ 4/.tm to 6/im spec- 
tral region (for furth er details, see iMarengo et al.l 120101 : 
iFreedman etaI1l2011l:[Scowcroft et al.ll2011l) . Infact. the 



^ For co nsistency, we only use the linear slopes, b^iu 
IBono et al] l(201fl,). 



from 



Theoretical IRAC Band P-L and P-C Relations 



3 



slight increase of 4.5/im and 5.8/im band P-L slopes was 
shown in all model sets presented in Table [1] 

4. COMPARISON TO THE EMPIRICAL P-L 
RELATIONS 

Slopes of the current empirical IR AC band P-L re- 
lations are summarized in Table 2 of iNgeow fc Kanbuii 
llblO). These include six sets of P-L relations in our 
Galaxy and Magellanic Clouds. Briefly, GALl and 
GAL2 arc P-L slopes derived from the "old" and "new" 
IRSB distances, respectively, and GAL3 are slopes de- 
rived from eight Cepheids that possess Hubble Space 
Telescope (HST) parallax measurements. Details con- 
cerning these Galac tic P-L slopes can be found in 
iMarengo etldl ()2010D . The LMCl and LMC2 are the 
empir ical L MC P-L slopes ta ken from iMadore et al.l 
(|2009D and INgeow et al.l (|2009D . res pectively, and the 
SMC P-L slopes were adopted from Ngeow & Kanbur" 
IJOIO). As discussed in .Ngeow & Kanbur (2010), these 
six sets of P-L slopes can be grouped to two groups char- 
acterized by steeper (~ —3.46) and shallower (~ —3.18) 
slopes respectively. Both of the steeper and shallower 
slopes can be predicted from using Lx = 4:'nB?B\{T), by 
assuming the behavi or of B\ (T) at long wavelengths for 
the IRAC bands (seelFreedman et al.ll2008l : iNeilson et"all 
[2?)Tnl: INgeow et al.ll20inl for more details). 

In addition to these six sets of empirical P-L 
relations obtained fr om random-phase observations, 
iScowcroft et al.l (|2011[ ) have recently published LMC 
P-L relations based on ^ 80 Cepheids that possess 
mean magnitudes in 3.6^m and 4.5/im bands. These 
LMC Cepheids have been observed multiple times us- 
ing Spitzer, with 24 evenly spaced data points per light 
curves, hence accurate mean magnitudes can be ob- 
tained. Their adopted P-L relations in 3.6/im and 4.5/im 
bands are denoted as LMC3 in this paper. 

These empirical P-L relations can be compared to the 
synthetic IRAC band P-L relations from the s e lected 
model sets given in Tables [D and H IBono et al.l (pOlOl ) 
have compared the synthetic BVIJK band P-L slopes of 
these selected model sets to their empirical counterparts 
for Galactic and Magellanic Clouds Cepheids, and found 
they are generally in agreement. When comparing the 
P-L intercepts, three different values of LMC and SMC 
distance moduli (/iLMC,SMc)l3 respectively, were adopted 
(see the right panels of Figures H] and [S|) . These distance 
moduli covered a wide range of available distance moduli 
in the literature for the Magellanic Clouds. 

4.1. Comparison of the Galactic P-L Relations 

A comparison of the synthetic and empirical Galactic 
P-L relations in IRAC bands are presented in Figure [3l 
The left panel of this figure shows that the synthetic P-L 
slopes from 12 + \og{0/H) = 8.58 {Z = 0.01, Y = 0.26) 
and 12 -f \og{0/H) = 8.89 {Z = 0.02, Y = 0.26) 
model sets are in marginal agreement with the GAL2 
and GAL3 P-L slopes, but not for the GALl P-L slopes 
derived from the "old" IRSB distances. However, the 
synthetic P-L intercepts from these two model sets agree 

^ Extinction i s ign o red as it is negligible in IRAC bands 
(IFreedman et al.l [20C^ : IFreedman & Madord [20Tol : INgeow et aP 
120091) . 



well with the GAL3 P-L intercepts (albeit the large er- 
ror bars) , but disagree with the other tw o empirical P-L 
intercepts based on the IRSB methods. IMarengo et al.l 
llOlO) have compared the synthetic P-L relations from 
12 + \og{0/H) = 8.90 (Z = 0.02, Y = 0.28) model set to 
the three sets of empirical Galactic P-L relations. They 
are also included in Figure [3] T he comparison shown 
in this figure echoes the finding in IMarengo et al.l ()2010f ) 
that this set of synthetic P-L relations agree well with 
their empirical counterparts derived from the /fS'T' paral- 
laxes (GAL3), but not for the other two sets of empirical 
Galactic P-L relations. 

Determination of the Galactic P-L relations is expected 
to improve in the near future. The Carnegie Hubble Pro- 
gram will observe and measure the distances to about 
39 Galactic Cepheids in 3.6/im and 4.5/im, with 15 of 
them expected to have paral lax measurements from Gaia 
llFreedman fc Madord l2010f l . The improved empirical 
Galactic P-L relations are expected to be able to dis- 
criminate the synthetic P-L relations that are best at 
describing the observed P-L relations. 

4.2. Comparison of the LMC P-L Relations 

In Figure lU synthetic P-L relations from three model 
sets were compared to the empirical LMC P-L relations. 
The synthetic P-L slopes from 12 -\- log{0/H) = 8.47 
{Z — 0.008, Y = 0.25) model set are in good agree- 
ment with LMC2 P-L slopes from Nge ow et al.l (p009). 
Note that the model set with Z = 0.008 and Y = 0.25 
is generally adopted as a representative metallicity for 
LMC Cepheids. On the other hand, the synthetic P- 
L slopes from the 12 -h log(0/i?) = 8.35 (Z = 0.006, 
Y = 0.25) model set also agree with LMCl P-L slopes 
(from iMadore et all l2009( ) in 3.6/im and 4.5/im bands, 
but not in the two longer wa velength bands. Sim- 
ilarly, LMC3 P-L slopes from IScowcroft et all (|201lH 
agreed with the synthetic P-L slopes from the same 
model set. For the P-L intercepts, the right panels of 
Figure |4] show that the synthetic P-L intercepts from 
12 -I- log{0/H) = 8.35 model set match with the em- 
pirical results from LMC2 if the LMC distance modulus 
ilJ-LMc) is adopted to be 18.50 mag. The LMC2 P-L in- 
tercepts at 5.8/im and 8.0/im bands also matched to the 
synthetic P-L intercepts from 12 -\- \og{0 / H) = 8.47 and 
12 + log(0/iJ) = 8.17 model sets if u^.mc - 18.6 mag ^ 

The P-L slopes adopted from IScowcroft et al.l ()2011l ) 
are based on Cepheids w ith period rang e of 1.0 < 
log(P) < 1.8. Table 3 of IScowcroft et alJ (pOUl l also 
listed the P-L relations using Cepheids with 0.8 < 
log(P) < 1.8, at which the slopes of these P-L rela- 
tions (-3.31 ± 0.05 and -3.22 ± 0.05 at 3.6/im and 
4.5/im bands, respectively) are in better agreement to 
the LMC2 empirical P-L slopes from ' Ngeow et al.l ()2009f ) 
and the synthetic P-L slopes from Z = 0.004 and 
Z = 0.008 (both with Y = 0.25) model sets given in 
tje previous section. The synthetic P-L slopes using re- 
stricted period ranges, either for 1.0 < log(P) < 1.8 or 
0.8 < log(P) < 1.8, ho wever, do not re p roduc e the steep 
slopes as adopted in IScowcroft et al.l (|2011| ). In con- 
trast, these P-L slopes are shallower than the synthetic 
P-L slopes from the full period range (0.4 < log[P] < 
2.0). For examples, 3.6/im band P-L slopes for Z = 
(0.004, 0.006, 0.008) model sets with 1.0 < log(P) < 1.8 
are -3.153 ± 0.019, -3.223 ± 0.022 and -3.134 ± 0.017, 
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Fig. 2. — Synth etic P-L slopes as a function of wavelength for selected models. Theoretical slopes in BVIJK bands were adopted from 
IBono et al] | |20TO| . their Table 2). 



TABLE 1 

Slope of the theoretical IRAC band P-L relations at various metallicities. 





z 




Y 


\og{Z/X) 


12 + log{0/H)^ 


[Fc/H]'' 


AY/AZ"" 


3.6/im 


4.5/im 


5.8/im 


S.O^lm 





0004 





24 


-3.28 


7.16 


-1.50 


25 


-3.511 ±0.004 


-3.506 ± 0.004 


-3.513 ± 0.004 


-3.519 ±0.004 





001 





24 


-2.87 


7.56 


-1.10 


10 


-3.344 ± 0.009 


-3.323 ± 0.009 


-3.336 ± 0.009 


-3.352 ±0.008 





002 





24 


-2.58 


7.86 


-0.80 


5 


-3.510 ±0.007 


-3.485 ± 0.007 


-3.498 ± 0.007 


-3.517 ±0.007 





004 





25 


-2.27 


8.17 


-0.49 


5 


-3.292 ±0.005 


-3.253 ± 0.006 


-3.269 ± 0.006 


-3.298 ±0.005 





006 





25 


-2.09 


8.35 


-0.31 


3.3 


-3.365 ±0.007 


-3.304 ± 0.008 


-3.326 ± 0.008 


-3.371 ±0.007 





008 





25 


-1.97 


8.47 


-0.18 


2.5 


-3.244 ± 0.005 


-3.187 ±0.006 


-3.206 ± 0.006 


-3.248 ±0.005 





01 





26 


-1.86 


8.58 


-0.08 


3 


-3.361 ±0.007 


-3.299 ± 0.009 


-3.318 ± 0.008 


-3.364 ±0.007 





02 





25 


-1.56 


8.88 


4-0.22 


1 


-3.311 ±0.005 


-3.239 ± 0.006 


-3.258 ± 0.006 


-3.312 ±0.005 





02 





26 


-1.56 


8.89 


-1-0.23 


1.5 


-3.369 ± 0.006 


-3.304 ± 0.007 


-3.322 ± 0.006 


-3.372 ±0.006 





02 





28 


-1.54 


8.90 


4-0.24 


2.5 


-3.13 ±0.01^^ 


-3.04±0.01<= 


-3.07±0.01<= 


-3.12 ± 0.01=^ 





02 





31 


-1.53 


8.92 


-fO.26 


4 


-3.271 ±0.005 


-3.191 ± 0.005 


-3.213 ± 0.005 


-3.272 ±0.005 





03 





275 


-1.36 


9.08 


-1-0.42 


1.5 


-3.245 ± 0.006 


-3.171 ± 0.007 


-3.191 ± 0.006 


-3.240 ±0.006 





03 





31 


-1.34 


9.10 


-1-0.44 


2.7 


-3.179 ±0.004 


-3.093 ± 0.005 


-3.118 ± 0.004 


-3.167 ±0.004 





03 





335 


-1.33 


9.12 


-fO.46 


3.5 


-3.297 ± 0.003 


-3.210 ±0.004 


-3.235 ± 0.004 


-3.285 ±0.004 





04 





25 


-1.25 


9.19 


-1-0.53 


0.5 


-3.343 ± 0.004 


-3.268 ± 0.005 


-3.289 ± 0.005 


-3.336 ± 0.004 





04 





29 


-1.22 


9.22 


+0.56 


1.5 


-3.182 ±0.005 


-3.104 ±0.005 


-3.125 ± 0.005 


-3.172 ±0.005 





04 





33 


-1.20 


9.25 


-1-0.59 


2.5 


-3.206 ±0.002 


-3.129 ± 0.003 


-3.150 ±0.002 


-3.195 ±0.002 



^ Calculated from an online tool, h ttp : //astro . wsu . edu/models/calc/MIX . html^ assuming 12 -t- lo§,{0/ H)q — 8.66. 

^ AY/AZ — (Y — Yp)/Z is the relative helium enrichment ratio, where Yp — 0.23 fsee lFiorentino et al.|[2002l . and reference therein). 
Taken from lMarengo et al.l I I2010I1 . 



respectively. Similarly, the 4.5/im band P-L slopes for 
Z = (0.004, 0.006, 0.008) model sets are -3.081 ±0.022, 
-3.125 ± 0.025 and -3.046 ± 0.021, respectively. This 
suggests that the adopted period r ange could affect the 
derived P-L relations. Interestinglv. INeilson et aTl ()2010f ) 
demonstrated that the empirical P-L slopes based on 
I Ngeow et alj ()2009() data can be steepen and co nsis- 
tent to the P-L slop es from iMadore et all (|2009[ ). or 
IScowcroft etHI (|20Tll ). if a period cut of log(P) = 1.05 
is applied. 

4.3. Comparison of the SMC P-L Relations 

Synthetic P-L relations from three model sets, 12 + 

log{0/H) = 7.86 {Z = 0.002, Y = 0.24), 12 + 

\og{0/H) = 8.17 (Z = 0.004, Y = 0.25) and 12 + 

\og{0/H) = 8.47 (Z = 0.008, Y = 0.25), were 

compared to the empirical SMC P-L relations from 



INgeow fc Kanbuil (|2010l) . As can be seen from Figure [5l 
the synthetic P-L relations from 12 + log(0/-ff) = 8.47 
model set agree with the empirical P-L relations, if the 
assumed SMC distance modulus is 19.10 mag. This re- 
sult is odd because Z — 0.004 is generally adopted as rep- 
resentative of SMC metallicity. Though there is a spread 
in the metallicity of SMC Cepheids from spectroscopic 
measurements, the mean value is closer to Z — 0.004 
than to Z — 0.008. This occurrence is due to the fact 
that current nonlinear pulsation models predict a signif- 
icant steepening of P-L slopes when changing the metal- 
hcity from Z = 0.008 to Z = 0.004 (with a significant 
reduction of the dependence at still smaller metal con- 
tents) whereas the empirical relations adopted in this 
paper tend to suggest almost the same slopes for LMC 
and SMC. For 12 + log{0/H) = 8.17 model set, the syn- 
thetic P-L slopes agree with the empirical P-L relations 
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TABLE 2 

Intercept of the theoretical IRAC band P-L relations at various metallicities. 





z 




Y 




12 + XoMIHY 


[Fc/HP 

L / J 




3.6/jm 


4.5//m 


5.8^m 


8.0/^m 





0004 





24 


-3.28 


7.16 


-1.50 


25 


-2.410 + 0.004 


-2.418 + 0.004 


-2.418 ± 0.004 


-2.424 + 0.004 





001 





24 


-2.87 


7.56 


-1.10 


10 


-2.527 + 0.010 


-2.547 + 0.011 


-2.542 ± 0.010 


-2.539 + 0.010 





002 





24 


-2.58 


7.86 


-0.80 


5 


-2.341 + 0.006 


-2.356 ± 0.006 


-2.353 ± 0.006 


-2.352 + 0.006 





004 





25 


-2.27 


8.17 


-0.49 


5 


-2.718 + 0.005 


-2.740 ± 0.006 


-2.737 ± 0.006 


-2.730 + 0.005 





006 





25 


-2.09 


8.35 


-0.31 


3.3 


-2.588 + 0.007 


-2.621 ± 0.008 


-2.615 ± 0.007 


-2.599 + 0.007 





008 





25 


-1.97 


8.47 


-0.18 


2.5 


-2.702 + 0.005 


-2.724 ± 0.006 


-2.723 ± 0.006 


-2.714 + 0.005 





01 





26 


-1.86 


8.58 


-0.08 


3 


-2.577 + 0.007 


-2.594 ± 0.008 


-2.595 ± 0.008 


-2.590 + 0.007 





02 





25 


-1.56 


8.88 


+0.22 


1 


-2.608 + 0.005 


-2.609 ± 0.005 


-2.617 + 0.005 


-2.621 + 0.005 





02 





26 


-1.56 


8.89 


+0.23 


1.5 


-2.596 + 0.005 


-2.598 ± 0.006 


-2.605 ± 0.006 


-2.607 + 0.005 





02 





28 


-1.54 


8.90 


+0.24 


2.5 


-2.66 + 0.01^^ 


-2.67 + 0.01<= 


-2.69 ± 0.01<= 


-2.68 ± 0.01=^ 





02 





31 


-1.53 


8.92 


+0.26 


4 


-2.608 + 0.004 


-2.627 + 0.005 


-2.629 ± 0.005 


-2.620 + 0.005 





03 





275 


-1.36 


9.08 


+0.42 


1.5 


-2.630 + 0.006 


-2.619 ± 0.006 


-2.629 ± 0.006 


-2.647 ± 0.006 





03 





31 


-1.34 


9.10 


+0.44 


2.7 


-2.637 + 0.004 


-2.639 ± 0.005 


-2.644 ± 0.004 


-2.660 ± 0.004 





03 





335 


-1.33 


9.12 


+0.46 


3.5 


-2.508 + 0.004 


-2.514 + 0.004 


-2.518 ± 0.004 


-2.532 + 0.004 





04 





25 


-1.25 


9.19 


+0.53 


0.5 


-2.545 + 0.004 


-2.525 ± 0.004 


-2.537 + 0.004 


-2.563 + 0.004 





04 





29 


-1.22 


9.22 


+0.56 


1.5 


-2.695 + 0.004 


-2.680 + 0.005 


-2.690 + 0.005 


-2.716 + 0.004 





04 





33 


-1.20 


9.25 


+0.59 


2.5 


-2.592 + 0.002 


-2.580 + 0.003 


-2.589 ± 0.002 


-2.615 + 0.002 



^ Calculated from an online tool, http://astro.wsu.edu/models/calc/MIX.htinl assuming 12 + log(0///)o — 8.66. 

^ Ay/AZ — (y — Yj)) / Z is t he relative helium enrichment ratio, where Yp — 0.23 fsee lFiorentino et al.|[2002l . and reference therein). 
Taken from lMarengo et al.l I I2010I) . 



<U -3.4 

a 

o 



- GALl 
o GALS 

□ GAL3 



Z = 0.01, Y = 0.26 

Z = O.Oa, Y=0.26 

Z = 0.Q2, ¥ = 0.28 



12 + log(0/H) = a.58 
12+log(0/H) = 6.89 
12 + log(0/H) = a.90 



Wavelength [micron] 



" GALl 
o GAL2 
□ GAL3 



Z = 0.01, Y=0.26 

Z=0.02, Y=0.26 

Z = 0.02, Y=0.28 



12 + log(0/H)=8.5B 
12+log(0/H)=8.B9 
12 + log(0/H)=8.90 




Wavelength [micron] 



Fig. 3. — Comparison of the empirical Galactic P-L relations, adopted from lMarengo et al.l 1)20101 ). to the selected synthetic P-L relations 
given in Table [T] (left panel) and [2] (right panel). Note that for better visuahzation, data points for GALl and GAL3 have been shifted 
shghtly in wavelength. See the text for the definition of GALl, GAL2 and GAL3. 



at 5.8/.tni and 8.0/im bands, though the 3.6/.tni and 4.5/im 
P-L slopes are slightly off from the empirical counter- 
parts. 

5. DEPENDENCY OF THEORETICAL P-L 
RELATIONS ON METALLICITY 

As mentioned in the Introduction, the IRAC band P- 
L relations are expected to be insensitive to metallicity. 
However, based on the residual analysis from multi-band 
empirical P-L relations for a number of LMC Ccpheids 
with [Fe/iJ] measurements, Freedman & Madorc (201l|) 
suggested the mid-infrared P-L relation could be mildly 
depending on metallicity. Therefore, possible metallicity 
dependence of the IRAC band P-L relations is investi- 
gated in this section using the synthetic P-L relations 
presented in Tables [T] and [H Figure |6] and [7] present the 



TABLE 3 

F-TEST RESULTS FOR SYNTHETIC P-L 
RELATIONS AS A FUNCTION OF 
METALLICITY. 



Band 


P-L Slope 


P-L Intercept 


3.6/^m 


15.3 


5.68 


4.5/^m 


24.3 


3.47 


5.8/im 


23.3 


4.34 


8.0Mm 


17.2 


6.46 



Note. — We only list the results for 
12 + \og{0/H). The _F-test results for 
[Fe/i?] are very similar and hence not 
listed in this Table. 
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Fig. 4. — Same as Figure [S] assuming three different values of the LMC distance modulus. These distance moduli roughly cover the 
available LMC distance modulus in literature. Note that for better visualization, data points for LMC3 have been shifted slightly in 
wavelength. 
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Fig. 5. — Same as Figure [S] assuming three different values of the SMC distance modulus. These distance moduli roughly cover the 
available SMC distance modulus in literature. 



synthetic P-L slopes (left panels) and intercepts (right 

panels) as a function of 12 + log(0/ijjf|, suggesting the 
IRAC band P-L relations may be sensitive to metallicity 
and/or helium abundance of the Cepheids. 

Well known statistical tests can be used to test the 
possible metallicity dependence of the synthetic P-L re- 
lations as shown in Figures |6] and [T) We considered a 
null hypothesis that the IRAC band P-L relations are 
independent of metallicity, which can be represented by 
a constant regression model in the form oi A ~ gq, where 

^ Plots for the synthetic P-L relations as a function of log{Z/X) 
and [Fe/H] are similar to these figures, hence omitted from the 
present paper. 



A stands for either the P-L slopes or intercepts. The al- 
ternate hypothesis is that there is a linear dependency of 
metallicity on P-L relations represented by a linear re- 
gression model, A = oq + oiB where B = 12 + \og{0/H) 
or B — [Fe/H]. Graphic representations of these re- 
gression models are presented in Figure |51 It is worth 
to point out that these regression models are used to 
test the dependency of metallicity on synthetic P-L rela- 
tions, and do not represent the actual linear dependency 
of metallicity on P-L relations, since the metallicity ef- 
fect, if present, will not be a simple linear relation (as 
evident in Figure [5] and [7]) . It is well-known in statistical 
literature (for example, see lKutner et al.l 12005') that the 
i^-test can be applied to test if additional parameter (in 
our case, the oi) is needed in the regression model. By 
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Fig. 6. — Synthetic P-L relations as a function of 12 -|- log(0/H), separated by the helium abundance (Y). The left and right panels are 
for the P-L slopes and intercepts, respectively. 
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Fig. 7. — Synthetic P-L relations as a function of 12 -|- log(0/H), separated according to the AY/AZ values. Note that the open symbols 
are for AY/AZ > 3, and filled symbols are for AY/AZ between 1 and 3 (inclusive). 



adopting a = 0.05, the null hypothesis can be rejected if 
F > 4.54. The F-test results are summarized in Tabled 
showing that the null hypothesis can be rejected in most 
cases, except for 4.5/ini and 5.8/im band P-L intercepts. 
These suggested that the synthetic P-L relations may not 
be independent of mctallicity. 

6. THE SYNTHETIC IRAC BAND P-C RELATIONS 

Pulsators from various model sets, as described in Sec- 
tion 2, can also be used to construct the synthetic P-C 
relations in the IRAC bands. For brevity, colors from 
3 . 6 fim and 4 . 5/im bands are denoted as [3 . 6] — [4 . 5] , and so 
on. Two examples of the synthetic P-C relations are pre- 
sented in Figure O Results of the synthetic P-C relations 
are summarized in Table |4l As in the case of the syn- 
thetic P-L relations, pulsators with 0.4 < log(P) < 2.0 



were used to fit the P-C relations. The P-C slopes and 
intercepts do not deviate by more than 0.006 if all the 
pulsators were included. The synthetic P-C slopes and 
intercepts were plotted as a function of 12 + \og{0/H) 
in Figures [TUl and [TTl respectively. From these figures, it 
is clear that P-C relations exist for certain combinations 
of IRAC band filters and mctallicity, and not all of the 
synthetic P-C relations are independent of metallicity. 
The P-C relations that are independent or insensitive to 
metallicity are the [3.6] - [8.0] and [4.5] - [5.8] P-C rela- 
tions, especially for those with 12 + \og{0/H) < 8.9. 

6.1. Comparison to Empirical P-C Relations 

iMarengo et all (|20Tol) found a significant [3.6] - [8.0] P- 
C relation for the Galactic Cepheids, with the expression 
of [3.6] - [8.0] = 0.039(±0.008)log(P) - 0.058(±0.014). 
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Fig. 8. — Fitting of the constant (upper panels) and linear (lower panels) regression models to the synthetic P-L slopes (left panels) and 
intercepts (right panels). The dashed lines represent the Icr boundary of the fitted regressions. The plots for [Fe/H] are very similar to 
this figure. 



This empirical non-zero P-C slope is in contradiction 
with the synthetic P-C slopes given in Table SI at which 
the synthetic [3.6] — [8.0] P-C slopes are close to zero. 
Disagreement was also found for the P-C intercepts. 
In Figure [121 P-C relations for 26 fundamental model 
Galactic Cepheids were compared to the synthetic P- 
C relations from selected model sets. IRAC band pho- 
tometry for these Ga lactic Cepheids were adopted from 
iMarengo etHI (|2010( ). This fi gure shows that only the 
[3.6] - [4.5] and [5.8] - [8.0] synthetic P-C relations barely 
agreed with the observed P-C relations. Disagreements 
between empirical and synthetic P-C relations may due 
to the small number of Cepheids in the sample, large pho- 
tometric errors, colors for all Cepheids are not be mea- 
sured at mean light, presence of circumstellar envelops 
especially for longer period Cepheids, uncertainties in 



the adopted model atmospheres or any combinations of 
these. 

In addition to Galactic Cepheids, iScowcroft et al.l 
()2011| ) also reported the finding of [3.6] - [4.5] P-C re- 
lation based on LMC Cepheids, within the period range 
of 1.0 < log(P) < 1.8, that possess multi-epoch observa- 
tions. The left panel of Figure [13] compares the P-C rela- 
tion for these Cepheids and the three synthetic P-C rela- 
tions given in TablelH where the [3.6] — [4.5] colors are ob- 
tained fro m the 3.6/xm and 4 .5/im band mean magnitudes 
taken from IScowcroft et al.l (!2011). It is clear that these 
synthetic P-C relations do not agree with the empirical 
P-C relation of [3.6 ] - [4.5] = -0.087(±0 012) log(P) + 
0.092(±0.013) from IScowcroft et al.l (|2QTTI ). However, as 
in Section 4.2, this empirical P-C relation should also be 
compared to the synthetic P-C relations based on pul- 
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TABLE 4 

Synthetic IRAC band P-C relations at various metallicities^. 



Z Y [3.6] - [4.5] [3.6] - [5.8] [3.6] - [8.0] [4.5] ~ [5.8] [4.5] - [8.0] [5.8] - [8.0] 



P-C Slopes 



0.0004 


0.24 


-0.005 


0.002 


0.008 


0.007 


0.012 


0.006 


0.001 


0.24 


—0.021 


—0.008 


0.008 


0.013 


0.029 


0.016 


0.002 


0.24 


-0.024 


-0.011 


0.008 


0.013 


0.032 


0.019 


0.004 


0.25 


-0.039 


-0.023 


0.006 


0.016 


0.045 


0.029 


0.006 


0.25 


-0.061 


-0.039 


0.006 


0.023 


0.067 


0.045 


0.008 


0.25 


-0.057 


-0.038 


0.004 


0.019 


0.061 


0.042 


0.01 


0.26 


-0.062 


-0.043 


0.003 


0.019 


0.065 


0.045 


0.02 


0.25 


-0.072 


-0.053 


0.001 


0.019 


0.073 


0.054 


0.02 


0.26 


-0.065 


-0.047 


0.003 


0.018 


0.067 


0.050 


0.02 


0.28 


-0.09*^ 


-0.06*^ 


-0.01*= 


0.03*^ 


0.08*^ 


0.05'' 


0.02 


0.31 


-0.079 


-0.058 


0.001 


0.021 


0.081 


0.059 


0.03 


0.275 


-0.074 


-0.054 


-0.005 


0.020 


0.069 


0.050 


0.03 


0.31 


-0.086 


-0.061 


-0.012 


0.025 


0.074 


0.049 


0.03 


0.335 


-0.088 


-0.063 


-0.012 


0.025 


0.075 


0.051 


0.04 


0.25 


-0.075 


-0.055 


-0.007 


0.020 


0.068 


0.048 


0.04 


0.29 


-0.078 


-0.057 


-0.010 


0.022 


0.068 


0.047 


0.04 


0.33 


-0.078 


-0.057 


-0.011 


0.021 


0.067 


0.045 


P-C Intercepts 


0.0004 


0.24 


0.008 


0.008 


0.014 


-0.000 


0.006 


0.006 


0.001 


0.24 


0.020 


0.015 


0.012 


-0.005 


-0.008 


-0.003 


0.002 


0.24 


0.015 


0.012 


0.011 


-0.003 


-0.004 


-0.001 


0.004 


0.25 


0.022 


0.019 


0.012 


-0.003 


-0.010 


-0.007 


0.006 


0.25 


0.034 


0.027 


0.011 


-0.006 


-0.023 


-0.016 


0.008 


0.25 


0.022 


0.021 


0.012 


-0.001 


-0.011 


-0.009 


0.01 


0.26 


0.017 


0.019 


0.013 


0.002 


-0.004 


-0.006 


0.02 


0.25 


0.001 


0.008 


0.013 


0.007 


0.012 


0.005 


0.02 


0.26 


0.002 


0.008 


0.011 


0.007 


0.010 


0.003 


0.02 


0.28 


0.01*^ 


0.03*= 


0.02*^ 


0.02*= 


O.Ol'' 


-0.01*= 


0.02 


0.31 


0.019 


0.022 


0.012 


0.002 


-0.007 


-0.009 


0.03 


0.275 


-0.011 


-0.001 


0.017 


0.010 


0.028 


0.018 


0.03 


0.31 


0.002 


0.006 


0.022 


0.004 


0.021 


0.016 


0.03 


0.335 


0.006 


0.010 


0.023 


0.004 


0.017 


0.013 


0.04 


0.25 


-0.020 


-0.009 


0.018 


0.011 


0.038 


0.027 


0.04 


0.29 


-0.015 


-0.006 


0.020 


0.010 


0.035 


0.026 


0.04 


0.33 


-0.012 


-0.003 


0.022 


0.009 


0.035 


0.026 



^ Errors of the P-C slopes and intercepts are ignored as they are generally less than 0.003. 
^ Derived from the P-L relations. 
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Fig. 9. — Two examples of the synthetic P-C relations in the 
IRAC band. Top and bottom panels show the example of a flat P- 
C relation and a P-C relation with largest slope, respectively. The 
corresponding pulsational model sets are given in the upper-left 
corners. 



TABLE 5 

Synthetic [3.6] — [4.5] P-C relation for 

THREE MODELS SETS WITH PULSATORS OF 

1.0 < log(P) < 1.8. 



z 


Y 


P-C Slope 


P-C Intercept 


0.004 


0.25 


-0.072 ± 0.003 


0.066 ± 0.004 


0.006 


0.25 


-0.098 ± 0.004 


0.082 ±0.005 


0.008 


0.25 


-0.088 ± 0.004 


0.065 ±0.005 



sators within the period range of 1.0 < log(P) < 1.8. 
These synthetic P-C relations for the three model sets 
are summarized in Table [5l and compared to the empir- 
ical P-C relation in the right panel of Figure [131 Table 
[5] shows that the synthetic P-C slope from (Z = 0.008, 

Y = 0.25) model set is in excellent agreement to the 
empirical P-C slope, although the synthetic P-C inter- 
cept is 2cr smaller than the empirical counterpart. On 
the other hand, synthetic P-C relation from (Z = 0.006, 

Y = 0.25) model set is also consistent with the empirical 
P-C relation. 

6.2. The Wesenheit Function 
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Fig. 10. — Slopes of the synthetic P-C relations as a function of 12 + log(0/-ff). The dashed lines represent y = in order to guide the 
eyes, and not the fitting to the data. 
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Fig. 11. — Same as Figure [Tol but for the intercepts of the synthetic P-C relations. 



Existence of the P-C relations suggested the Wesen- 
heit function can be formulated in the IRAC bands. In 
optical bands, dispersion of the Wesenheit function is 

2 to '-^ 3 times s maller than the opti c al P-L relations 
(for e xample, see iFouaue et al.l l2007t iSoszvnski et al.l 
I2008t iNgeow et all 120091 and reference therein). The 
dispersions of the synthetic Wesenheit function, in the 
form ofW = m3.6Mm - 3.40 x ([3.6] - [4.5]), from aU 
model sets are on average about ~ 1.5 larger than 
the dispersions from the synthetic 3.6^m band P-L re- 
lations. Also, it is not worth to derive the IRAC 
band Wesenheit function because Wesenheit function is 
reddening-free by definition. On the other hand, extinc- 
tion c an almost be ignored for the IRAC band P-L rela- 
tions (iFreedman et a.1.1 l2008HFreedman fc Madord 120101: 
INgeow et al.ll2009l : IMarengo et al.ll2010l: IPreedman et al.l 



120111 ). Therefore, there is no net gain by using the We- 
senheit function in the IRAC bands. 

7. DISCUSSION AND CONCLUSION 

Synthetic IRAC band P-L relations based on a se- 
ries of pulsation models were investigated in this paper. 
These synthetic P-L relations were compared to their 
empirical counterparts, and the possible metallicity de- 
pendency of the IRAC band P-L relations was exam- 
ined. For the former part, selected sets of synthetic IRAC 
band P-L relations show agreement to the empirical P- 
L relations derived from Galactic and Magellanic Cloud 
Cepheids. The BVIJK synthetic P-L relations for these 
selected model sets also agreed with their optical and 
near-infrared counterparts, as presented in i Bono et al.l 
(|201Q) . for the Galactic and Magellanic Clouds P-L rela- 
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Fig. 12. — Comparison of the P-C relations for Galactic Ce plieids and the synthet ic P-C relations from three selected model sets given 
in Table |3] Photometry for Galactic Cepheids is adopted from IMarengo et al.l 120101) . 




Fig. 13. — Comparison of the P-C relation for LMC Cepheids based on the data presented in lScowcroft et at] 1)20111 ) and the synthetic 
P-C relations from three selected model sets. Left and right panels show the synthetic P-C relations using pulsators in the period range of 
0.4 < log(P) < 2.0 and 1.0 < log(P) < 1.8, respectively. 
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Fig. 14. — Comparison of the P-L intercepts between the em- 
pirical P-L relations and the synthetic P-L relations from selected 
model sets. The top panel shows the comparison for LMC, with 
synthetic P-L relations selected from {Z = 0.008, Y = 0.25) model 
set. The bottom panel shows the comparison for SMC, with syn- 
thetic P-L relations selected from {Z = 0.004, Y = 0.25) model 
set. 

tions. For the metallicity dependency of IRAC band P-L 
relations, plots for all of the synthetic P-L relations as a 
function of metallicity revealed that the IRAC band P-L 
relations may not be independent of metallicity. This re- 
sult is also supported by statistical i^-test. On the other 
hand, current empirical IRAC band P-L relations based 
on three galaxies, either the P- L relations in "steep" o r 
"shallow" groups (see Table 2 of lNgeow fc Kanburi[2mol ). 
suggested the IRAC band P-L relations may not be de- 
pending on metallicity, or have a weak dependency on 
metallicity. It is clear that more empirical determina- 
tions of the IRAC band P-L relations are needed, espe- 
cially for galaxies with low metallicity, from future ob- 
servations with JWST. 

Figure [HI shows the intercepts for LMC and SMC em- 
pirical P-L relations relative to the two adopted syn- 
thetic P-L relations, which is equivalent to derive the 



distance moduli to the Magellanic Clouds using the syn- 
thetic P-L relations. The resulted distance moduli for 
LMC and SMC are higher than the values commonly 
adopted in literature (e.g., 18.5 mag and 19.0 mag for 
LMC and SMC, respectively). This is due to the al- 
ready mentioned limitation of the adoption of a canoni- 
cal M-L relation. However, previous theoretical compu- 
tations of current nonlinear conyectiy e models (see, e.g. , 
iBono et al.l [200l ICamito et al.l [200l iBono et all 12001 
and references therein) have shown that by relying on 
non-canonical models based on an M-L relation brighter 
than the canonical one by 0.25 dcx the inferred distance 
moduli are shorter than the values obtained in the canon- 
ical scenario by 0.15-0.2 mag depending on the filters. 
This implies that adopting non-canonical relations in the 
comparison with empirical Spitzer data we would have 
obtained shorter distance moduli by 0.15-0.2 mag for 
each selected chemical composition, in better agreement 
with most recent and adopted values in the literature. 

Finally, synthetic P-C relations in the IRAC band were 
also derived and compared to the Galactic and LMC 
empirical counterparts. In general, disagreements were 
found between the synthetic and empirical P-C relations. 
However, the synthetic P-C relations are in agreement 
with the empirical LMC [3.6] - [4.5] P-C relation if a 
period range of 1.0 < log(P) < 1.8 is adopted when con- 
structing the synthetic P-C relations. Observations of a 
large number of Cepheids with JWST may help in re- 
solving the discrepancy of the IRAC band P-C relations. 
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